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Abstract. An overview of the flavour problem is presented, with emphasis on the theoretical 
efforts to find a satisfactory description of fermion masses and mixing angles. 


1. Introduction 

The origin of the parameters in the flavour sector of the Standard Model (SM), minimally 
extended to include massive neutrinos, is one of the most enigmatic questions in particle physics. 
Out of the 22 (20 if B-L is conserved) independent low-energy parameters Ti, which with some 
abuse of language can be called Yukawa couplings, 18 have been measured. Of the remaining 
four parameters, the absolute scale of neutrino masses is constrained in a limited range, the 
leptonic Dirac CP-violating phase starts to be constrained by global fits while the two possible 
Majorana phases are still unknown. A considerable effort has been devoted to search for a more 
economic description, perhaps related to a new principle, such as the gauge principle. Gauge 
invariance and renormalizability allow to describe strong and electroweak interactions of three 
copies of fifteen different fermion species in terms of only three parameters. Nothing similar 
exists so far in the flavour sector and we usually refer to this as the flavour puzzle. Another 
aspect of flavour is related to the new particle threshold around the TeV scale predicted by 
all SM extensions addressing the hierarchy problem. Once new TeV particles transforming 
non-trivially in flavour space are introduced, it is very difficult to maintain the almost perfect 
agreement between predictions and observations that reigns in the SM. New sources of flavour¬ 
changing neutral currents (FCNC) and CP violations appear and the task is to keep them at 
an acceptable level. This is what we commonly mean by flavour problem, to distinguish it from 
the flavour puzzle. In this short review the focus will be on the hrst aspect, I will comment only 
shortly on the second one. Also, I’m not aiming at reviewing all the existing models, but rather 
at revisiting some main ideas, guided by my own prejudices. Several aspects that have been left 
aside or only briefly mentioned in this paper are covered in a number of reviews PEIEIII!. In 
particular I will not do justice to the vast literature of fermion masses in grand unified theory 
mEiE], nor to that discussing the flavour puzzle in the framework of string theory 019]. 

There are different approaches to the flavour puzzle, with many intermediate possibilities. 
We may take a reductionist perspective: the Yukawa couplings Vi should be deduced from 
first principles. We postulate the existence of a fundamental theory from which Vi can be 
uniquely determined. Either by proceeding directly from the candidate theory or by appealing 
to some symmetry or dynamical principle, Vi are then computed in terms of a small set of input 


parameters. Probably the most striking fact about this program is that nothing approaching a 
standard theory of exists, despite the decades of experimental progress and theoretical efforts. 
In another approach a major role is played by chance. There are many variants and practical 
implementations of this strategy. The Yukawa couplings Ti are typically mapped to a large 
number of order-one parameters that are considered as irreducible unknowns, like in models 
with Froggatt-Nielsen abelian flavour symmetries or with fermions living in extra dimensions. 
Also the simplest version of partial compositeness falls into this class. By scanning the order-one 
parameters we get probability distributions for masses and mixing angles. Alternatively we start 
from a fundamental theory, like string theory, which possesses a vast landscape of solutions, with 
no privileged ground state. The observed Yukawa couplings become environmental quantities 
and cannot be predicted, like the relative sizes of the solar planetary orbits m- We are allowed 
to ask much less ambitious questions. For instance, if we have knowledge of the statistical 
distribution of Ti in an hypothetical multiverse where the laws of physics follow our fundamental 
theory, we can ask how typical are the Yukawa couplings that we observe. Conversely, barring 
anthropic selections, we might assume that the observed Ti are typical and try to deduce 
information on the statistical distribution of Ti in the multiverse [111 US]. Such a variety of 
open possibilities shows how far we are from the solution of the puzzle and even from identifying 
the most relevant questions to be addressed. 

Most of the parameters Ti are dimensionless and in a dynamical theory of flavour we have 
essentially no clue about the characteristic scale Af. If active neutrinos are Majorana particles 
and B-L gets violated at a scale A, then « 0.05 eV strongly suggests a very large 

A. However in general A and Af are independent from each other. Thus there is no clear 
relation between Af and other possible particle physics thresholds such as the TeV scale, relevant 
to the gauge hierarchy problem, or the grand unified scale. This makes it more difficult to 
identify unambiguous signatures to confirm or rule out a given model of fermion masses and 
mixing angles. For instance, the extrapolation of the Yukawa couplings from the scale Af down 
to low-energies where they are measured can involve new particle threshold and/or unknown 
parameters, thus affecting our ability to test the high-energy theory. 

2. Lessons from the quark sector 

A first useful observation is that ratios of charged fermion masses and quark mixing angles can 
be represented by powers of the Cabibbo angle. Using A = 0.22 we have 
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where all masses have been renormalised at the scale mz- It is well-known that also the elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix can be expressed in terms of powers 
of A: 

|Kd|~l \Vus\-\ lUefel^A^ iKbl-A^-A^ . (3) 

For comparison, in the lepton sector, where the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
mixing matrix is denoted by U, we have all \Ufi\ of order one, except for {Uesl which is of 
order A. The ratio between the solar and the atmospheric neutrino squared mass differences 
Am 2 i/\Am‘^^\ is of order A^. Focussing on the quark sector, in a pioneering work [T3] Froggatt 
and Nielsen observed that all the small dimensionless parameters of the quark sector such as 
the quark mass ratios and the CKM mixing angles can be interpreted as powers of the breaking 
parameter of a flavour symmetry. In this case the flavour symmetry group Gf is abelian, 
Gf = U{l)F]sf. A scalar field cp, carrying by convention a negative unit of the abelian charge 



FN, develops a vacuum expectation value (VEV) that can be parametrized as 


X = {cp)/Af<l FN{ip) = -l . (4) 

Quarks carry non-negative U(l)^^ charges (the case with charges of both signs can be discussed 
as well) 

FN{Xi)>0 (X, = . (5) 

Under these assumptions the quark Yukawa couplings yu^d are given by: 

yu = FucYuFq , yd = FdcYdFq , ( 6 ) 

where Yu^d are complex matrices with entries of order one, undetermined by the U(l)^^ 
symmetry, while Fx are real diagonal matrices, completely specified in terms of A by the charges 
FNiXi): 
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The small quark mass ratios and quark mixing angles originate from the hierarchical structure 
of the matrices Fx- Indeed, by taking FN{qi) > FN{q 2 ) > FN{q 3 ) > 0 we get 

iVu,dh « ^ < 1 ii<j) ( 8 ) 

for the matrices Vu 4 defining the CKM mixing matrix Vckm = VuVd- Independently from the 
specific charge choice, this framework predicts 

Vud-Vcs^Vtb-0{l) Vub-Vtd^VusXVa, , (9) 

the last equality being correct within a factor of two. With A ~ 0.2, the correct order 
of magnitudes of the Vckm matrix elements can be reproduced by choosing, for instance, 
FN{q) = (3, 2,0). The correct order of magnitudes of the quark mass ratios can be reproduced 
by choosing, for example 

FiV(g) = (3,2,0) FfV(nQ = (4,2,0) FN{(f) = {1 + r,r,r) , (10) 

r being a non-negative integer. If there is only one Higgs doublet, then we need r close to 
2 to match the ratio mtimb- If two Higgs doublets are present, other choices are possible by 
varying tan/3 = Vujvd- Several aspects of this class of models have been discussed in refs. 

[nmaoiiiiTiiiBiiiHiEoiEi]. 

The construction relies on a spontaneously broken abelian flavour symmetry, but the final 
results (|6]|^ are valid in a more general context, where no symmetry is present to start with. 
A simple example is provided by a model with an extra spatial dimension, compactified on an 
orbifold jZ^ to allow for 4D chiral fermions. The Lagrangian for a 5D spinor '^{x,y) reads: 

C = + 

= -h -F ... (11) 

where the mass m should be odd under the Z 2 parity sending y into —y. A possible choice is 

m = Me{y) , (12) 


M being a real constant and e{y) the periodic sign function. The 5D spinor has left (L) and 
right (R) chiralities in four dimensions 
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(13) 


with opposite Z 2 parities, such that only the even component developes a massless (zero) mode. 
Choosing, for instance, even and odd, the equation satisfied by the zero mode of is: 

a,T0+Me(2/)'1/0=0 . (14) 


The solution has an exponential dependence on y 
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where the first factor provides the correct normalization. The zero mode is localized near 
y = 0{TrR) for M > 0(< 0). In the limit M = 0 the zero mode becomes flat in y. A formally 
identical solution holds for the zero mode of /hi?, if we choose T l odd and even and we start 
from a 5D mass term with the opposite sign. If the Higgs field is strictly localized at one of 
the two branes, for instance the one at y = 0, the Yukawa interactions will be proportional to 
a Dirac delta 6{y) and we can reproduce the same pattern of Yukawa coupling of eq. (§ with 
matrices Fx now given by j22| 


Fx,= 



(16) 


where m and p are specified in terms of the bulk quark masses and the geometry of the extra 
dimension, see table The suppression factors Fx^ represent the values that the profiles of 


Table 1. Parameters pi and p in models with an extra dimension compactified on an interval. 
In a flat (warped) metric the Higgs field is localized on the brane y = 0 {y = R') and the 
ultraviolet cut-off is denoted by A {1/R). The fermions are described by five-dimensional spinors, 
with bulk masses Mi. In the warped case, when the framework is applied to the gauge hierarchy 
problem [23] , R, R! are length scales of the order of the inverse Planck mass and the inverse TeV 
scale, respectively. 
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the fermion zero modes take at the brane where the Higgs field is localized. On that brane, 
generic 0(1) Yukawa couplings with the bulk quark fields are postulated. The role of the 
Froggatt-Nielsen charges is here played by pi and p that determine the profiles along the extra 
dimension of the zero-mode wave functions: 




1 — 


Vii > 1 

1 < 1 

^/=|:eW2 <c-l 


y/pFx, 


{^i = 2p,p) . (17) 











There is no flavour symmetry: the hierarchical structure of quark masses and mixing angles is 
dictated by geometry in the compact space. 

Similarly, in the partial compositeness scenario [21], light fermions get hierarchical masses 
from the mixing between an elementary sector and a composite one. As a toy realization of 
this idea, consider a model where the composite sector contains, for each SM fermion, a pair of 
heavy fermions allowing a Dirac mass term of the order of the compositeness scale and a mixing 
term with the SM fields [25l [26] 

Cy = -u^AuU-d^AuD-Q^Agq 

-U^MuU - D^MdD - Q^MqQ 

-C/'=y4l>tQ) - D'^Ydi^^Q) - - {Q^^)YdD + h.c. (18) 

The first line represents the mixing between elementary and composite sector, the second line 
displays Dirac mass terms for the fermions of the composite sector and the third line shows 
the Yukawa interactions that, by assumption, are restricted to the composite sector alone and 
described by strong couplings Yu^d,Yu,d > 1- By integrating out the composite sector under the 
assumption Mi 3> u, we get low-energy Yukawa interactions for the elementary sector whose 
leading order (LO) terms have the structure given in eq. Q with matrices Fx parametrizing 
the elementary-composite mixing: 

Fu^ = AuM-^ , Fdc = AdMj^ , Fg = M-^Ag . (19) 

The same pattern arises when matter chiral multiplets Yj of the MSSM are coupled to a 
superconformal sector in some finite energy range [SZlEHlEe], from an ultraviolet (UV) scale A 
down to a lower scale Ac. Generic 0(1) Yukawa couplings Yij at the scale A 

w = XiYijXjH -|- .... ( 20 ) 

can undergo a sizable renormalization induced by the corrections to the Kahler potential. In the 
superconformal window the chiral multiplet Xi can have a large positive anomalous dimension 
7 j and the Kahler potential at the scale Ac becomes 

K = Y,Zi{Ac)x\Xi + ... (21) 

i 

where 

Z,{A,) = Z,{X) Zi{X)^l ( 22 ) 

Moving to a basis of canonical kinetic terms, the Yukawa couplings at the scale Ac are 
renormalized 


/A \ ^ 

Yg{A,) = FxY^gFx, 

and we find again the same pattern of eq. (©, without imposing any symmetry. 

In the previous examples the anarchical pattern of Y^^d may result in strong bounds on 
the scale of new physics Axp associated to particles carrying flavour quantum numbers and 
representing new sources of FCNC and/or CP violation. In the absence of a concrete realisation, 
it is difficult to estimate reliably the corresponding effects, also because in general the scale of 
new physics Ajvp and the scale of flavour physics Aj are independent from each other. A 
possibility is offered by a spurion analysis m, analogous to that prescribed by Minimal Flavour 


Violation (MFV)[3T]. To this purpose we assume that the new degrees of freedom have non¬ 
trivial flavour properties and that the flavour-violating effects are completely specified by the 
same spurions that are responsible for fermion masses and mixing angles. Moreover we assume 
that the dominant flavour-violating contributions admit an expansion in power series of the 
spurion fields. We start by noticing that the pattern of eq. Q is compatible with the flavour 
symmetry Gf = SU(3)^ x SU(3 )h^ with quarks transforming only under SU(3)^ as 

g = (3,1,1) = (1,3,1) d= = (l,l,3) . (24) 

The full symmetry Gf is explicitly broken by both the matrices Yu^d and Fx ■ However it can be 
formally restored by treating and Fx as non-dynamical spurion helds possessing suitable 
transformation properties. To this aim the Yukawa couplings should transform only under the 
“hidden” group SU(3 )h^: 

Y^ = {3,3,1)h Yrf = (3,l,3)H . (25) 

The suppression matrices Fx are the interface between SU(3 )h^ and SU(3)^, and they are 
assigned appropriate transformations under both factors to guarantee the invariance of the 
Yukawa interactions described by eq. ([^ under SU(3)^ x SU(3 )h^. The starting point of the 
spurion analysis is similar to that of MFV. Indeed the maximal flavour symmetry felt by quarks 
is SU(3)^, as in MFV. However there are more spurions than in MFV, the irreducible ones 
including now Fq, Fu^, Fd<:, Yu and Yd- One of the most dangerous effects originates from the 
effective operator 

‘ (5#)(sV) +... (26) 

^^NPYd) ^ 

(Y^) representing an average 0(1) coupling. The contribution of this operator to the CP- 
violating ex parameter is enhanced at the level of both the hadronic matrix element and the 
QCD corrections and sets one of the most stringent bounds on the scale of new physics Ajvp, 
see Table Assuming a generic 0(1) phase for the overall coefficient we need 
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(birsL)'^ 

1.4 X 10^ 

2.5 X 10^ 

5.0 X 10"^ 

1.7 X lO-'^ 

(ftp SL){bLSR) 
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8.8 X 10-® 

2.9 X 10"® 


Table 2. Bounds on dimension-six AF = 2 operators, from ref. |32[ I33| . The overall coefficient 
of the operators is cxp/Aj^p- The operator in eq. (26) is the one in the second row. 


{Yd) AxP > 20 TeV (27) 

not to spoil the SM prediction for ex- This, together with other constraints, suggests that a 
fully anarchical pattern in Yu^d is probably not tenable if new flavoured physics is present at the 
TeV scale [Hi] . 












When such a spurion analysis is applicable, the estimate of eq. (26) represents a sort of 
lower bound on the size of the expected effect and larger contributions are possible [3^. For 
example in supersymmetric extensions of the SM with a U(l)i;’Ar flavour symmetry the operator 
considered in eq. (26) receives contributions from box diagrams with squarks/gluino exchange 
that are typically larger than the one quoted in eq. (26). The reason is that in the \J{1)fn 
case the true flavour symmetry is much weaker than SU(3 )h^ x SU(3)^ and it allows sizable 
off-diagonal terms in both LL and RR blocks for the first two generations of the down squark 
mass matrix. For instance, with the charge assignment of eq. (10), the mass insertions { 5 i 2 )ll 
and { 812 )RR are both proportional to A and the operator {sd^){s‘^d) has an overall parametric 
suppression x A^/A^p, milder than the one in eq. (26). 


3. Prom quarks to leptons 

In the lepton sector we have no evidence for strong hierarchies in mixing angles or in neutrino 
masses. Hierarchy shows up at the level of charged lepton masses. In terms of the suppression 
factors Fx^ this means 

Pel < < Peg and Pq PS Fi^ k. Fi^ . (28) 

For example an acceptable set of charges is 

FAr(e=) = (4,2,0) FN{1) = {s+ t,s,s) (s>0,t = 0,l) . (29) 

Here we focus on Majorana neutrinos. In the context of a type I see-saw mechanism right- 
handed neutrinos have their own suppression matrices P^/c. Yukawa couplings and the 
mass matrix M of read 


= FycY^Fi 


ye = FecYeFi 


M — FycYcFycMQ 


(30) 


where Yi/,e,c are complex matrices with unknown entries of order one and Mq is a mass parameter. 
At low-energy the active neutrino mass matrix rriu is given by 

m, = -Fi (Yjy-^Y.) Fi uVMo , (31) 


with no dependence on the suppression matrices Fyc. 

A drastic realization of this picture is the framework of Anarchy [Ml EZl [Ml Ea SO], which 
corresponds to the case 

Po = P«2 = P/3 or ^ = 0 (32) 

In the anarchic framework the mass matrix for light neutrinos is 


/ 0(1) 0(1) 0(1) \ 

nil, = 0(1) 0(1) 0(1) mo 

V 0(1) 0(1) 0(1) ) 



(33) 


with undetermined order-one matrix elements. This implies mixing angles and nentrino mass 
ratios of 0(1), in rough agreement with the data. No special values for these quantities is 
expected. Indeed, before we knew Oiz from the experiments. Anarchy successfully anticipated 
values close to the upper bound at the time. Global fits of present data hint at deviations of 
the atmospheric mixing angle 023 from 7r/4. Today these indications are still weak, as shown by 
the instability of the best fit value against different fitting procedures. The persistence of these 
deviations in future tests would further strengthen the case for Anarchy. Anarchy represents 
an extreme possibility and milder realization of the relations (32) are possible. For instance, in 








the context of SU(5) grand unified models, with a Froggatt-Nielsen U(l)^^ abelian symmetry, 
neutrino masses and mixing angles can be reproduced, at the level of order of magnitudes, by 
several choices of the FN charges for the 5 multiplets hosting the lepton doublets, as shown in 
table FN charges for fermions in the 10 representations can be suitably chosen so that, by 
varying the unknown order-one parameters, reasonable distributions for charged lepton mass 
ratios, quark mass ratios and quark mixing angles are obtained [m im sa El ss]. A naive 
comparison of the distributions for neutrino masses and mixing angles with data do not appear 
to favor Anarchy over the other possible charge assignments, as can be seen from fig. I would 
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Figure 1. Probability distributions of r = sin @ 13 , tan^ 0i2, tan^ 023, from ref. 

|43j . within type I see-saw. The modulus (argument) of the complex random coefficients has 
been generated in the interval [0.5,2] ([0,27r]) with a flat distribution. For A and A = 0.2 
has been used, for H and A = 0.4 is taken. The shaded vertical band emphasizes the 

experimental 2cj window. 


Table 3. Possible choices of TA charges for the 5 representation in a class of SU(5) grand 
unified models, from ref. jl3] . The second column shows the value of the FN symmetry breaking 
parameter optimizing the fit to fermion masses and mixing angles. 
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personally find more appropriate to use the term Anarchy to denote the approaches giving rise 



























to the results ([^ and (30 31) where the absence of any special pattern resides in the matrices 
Y, rather than to indicate the special case defined in eq. (32). 


If this framework also comprises new flavoured particles at the TeV scale, severe bounds from 
lepton flavour violation (LFV) apply, under assumptions analogous to those spelled for the quark 
sector. The irreducible sources of flavour violation in the lepton sector include the matrices Y^, 
Fgc and Fi and LFV can occur even in the limit of vanishing neutrino masses. Notice that, 
though MFV cannot be extended unambiguously to the lepton sector [36], it predicts no LFV 
if neutrinos are massless since in this limit the only relevant spurion in the lepton sector is Y^, 
which can always be chosen diagonal. The dipole operator contributing to LFV is 


A2 




flU 


F^^''{Fe<^YeY}YeFi)Hn 


(34) 


The charged lepton mass matrix is proportional to {Fe^YeFi). In general the combinations 
{FgcYeFi) and (FecYeYeYeFi) are not diagonal in the same basis, not even in the case of universal 
Fi of eq. (32), and radiative decays of muon and tau are expected. Agreement with the most 

ey) < 5.7 x 10“^^, requires A^rp well above 10 TeV 


constraining upper bound, BR{^ 

As in the quark sector, a completely anarchical matrix Yf> and fiavoured physics at the TeV scale 
are difficult to reconcile. A sufficient condition for the absence of LFV is that Y)., Tgc and Fi are 
diagonal in the same basis, as suggested in some models. Alternatively we can look for special 
forms of these matrices |49l I50| . possibly dictated by some symmetry requirements. 


4. Realizations in grand unified theories 

A welcome feature of the above description is that it can be adapted to grand unified theories 
(GUT) where quarks and leptons are hosted in the same multiplet of the gauge group. In SU(5) 
the gauge symmetry requires 


Fn — Fl.c — Fpc — F-\ 


10 


Fi = Fdc = F^ 


Fyc — Fi 


(35) 


Our previous results, eqs (10) and (29), come very close to this requirement if we choose r = s 


and t = 1. If we accept a couple of tunings in the unknown 0(1) parameters Yu^d, we can 
force the equality (35) and still have a decent description of both the quark and lepton mass 


spectrum. As we have seen F^c drops from the low-energy quantities. It is instructive to 
consider also the ansatz iV oc 1. In this case the hierarchy among fermion masses is entirely 
due to Fiq. From eqs. ( 6|30 ) we see that mass ratios in the up-quark sector are the square of 
the respective mass ratios in the down-quark and in the charged lepton sectors, which is correct 
in first approximation. The large lepton mixing corresponds to a large mixing among F quarks 
ED, unobservable in SM weak interactions, but with possible observable effects if transferred 
from quarks to squarks in SUSY extensions of the SM [52|. A minimal model with Higgs bosons 
in the 5 representation would lead to the unrealistic relation i/g = yjj but the contributions 
from other Higgs representations or from non-renormalizable operators can solve this problem 
PUSH without altering the picture. 

At first sight this description does not seem to be compatible with an SO(IO) GUT. The most 
general renormalizable Yukawa interaction of three copies of fermion generations transforming 
as 16 of SO(10) reads 


£y = — 16j 


Y^hOn + Y:^.120 h + Y:^,.126h 




16,' -|- h.c. 


(36) 


The pattern of Yukawa couplings in eq. ([^ can also be thought to arise from a rescaling of the 
fermions fields, with the constraint that fermions belonging to a given irreducible representation 
of the gauge group have to undergo the same renormalization. By assuming that the matrices 















^ 10 ) ^ 120 ) ^126 have complex elements of order one and that the fields 16 undergo a wave function 
renormalization 


16 ^ Fi6l6 


(37) 


we see that all members of a 16 representation are affected in the same way. Even accounting for 
the Clebsch-Gordan coefficients arising from eq. (36) and the different overall scales associated 
with the Higgs VEVs we cannot reproduce the observed hierarchies of u, d and e masses. 

Such a discouraging starting point has been successfully modified in a construction by Kitano 
and Li [55], recently revisited in ref. [56] . 

The model is a SUSY SO(IO) GUT realized in a flat five-dimensional space time, the fifth 
dimension being compactified on an interval [0, ttR] whose inverse size is of the order of the 
GUT scale. The N = 1 5D SUSY corresponds to an Y = 2 4D SUSY, which is broken down 
to Y = 1 as a result of appropriate boundary conditions. The model comprises a 5D vector 
supermultiplet decomposing as a 4D vector 45v multiplet plus a 4D chiral 45$ multiplet. In the 
bulk there are also three copies of 5D hypermultiplets, equivalent to 4D chiral multiplets 16 and 
16^, with bulk masses Mi [i = 1,2,3). The boundary conditions allow zero modes only for 45v 
and 16. A Yukawa superpotential analogous to eq. (36) is localized at the brane y = 0. Prior 
to the SO(10) symmetry breaking, the wave functions of fermion zero modes evaluated at y = 0 
effectively drive a rescaling of the Yukawa couplings, as described by (37). The suppression 
factors in Eig are (see eq. 


•^ 16 i 



(38) 


where /r* = Mi/A and p = AttR. At this stage the Yukawa interactions are not yet able to 
reproduce the known hierarchies of u, d and e masses. The key ingredient of the model resides 
in the gauge interaction of the hypermultiplets. The 5D SUSY gauge interaction contains a 4D 
Yukawa interaction among 16j, 16*^^ and 45$, controlled by the gauge coupling yg, that can be 
combined with the bulk mass term: 


- 16? 


Mi - \/2y545$ 


16. 


(39) 


Q /9 

The chiral multiplet 45$ has no zero mode but can acquire a non-vanishing VEV, (45$) = , 

that breaks SO(IO) down to SU(5)xU(l)x- The Yukawa interaction of eq. (39) gives rise to 
effective SO(10)-breaking bulk masses: 


p:: = Pi- 


k = \/2y5u|'^^/A 


(40) 


where is the U(l)x charge of the different SU(5) components inside the 16 multiplet: 
Qx = (—1,-|-3, —5) for r = (10,5,1). We are back to the SU(5) case, see eq. (35): 


F -W 

1 / 1 _ e-2MrP 


r = (10,5,l) 


(41) 


but now the profiles Fr^ only depend on four free parameters: pi and k. Neutrinos are described 
within a type I see-saw mechanism, as in eq. (31), with masses for heavy Majorana neutrinos 


originating from the VEV of the SU(5) singlet in the 126// representation. 

The model contains many parameters of order one. After rephasing of the relevant fields 
there are 27 real parameters coming from the matrices Yio, Y 120 , Y 126 and 8 real parameters 
describing the embedding of the two light Higgs doublets within 10// and 120//. Despite the 
large number of parameters the agreement with data is not a priori guaranteed, since there are 
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Figure 2. The distributions of minimized for NO and 10 in neutrino masses and for 

tan/S = 50, from ref. [56]. 
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Figure 3. The predictions for varions observables obtained for Xmin/^ ^ fo of NO 
and tan/3 = 50, from ref. j56| . 


only 4 profile parameters to describe hierarchical mass ratios and mixing angles. Indeed a fit to 
an idealized set of 17 observables leads to a good agreement only for large values of tan^S, for 
both normal (NO) and inverted (10) neutrino mass ordering. 

A closer inspection reveals that fitting fermion masses and mixing angles in the 10 case 
requires a fine-tuning of the Yukawa parameters. By generating a large sample of random 
order-one Ynkawa parameters, the fit can be repeated by keeping at each iteration only 12 free 
parameters, 4 for the profiles and 8 for the relevant Higgs combinations. The distributions of the 
minimum over the number of degrees of freedom are shown in fig. [^for NO and 10. We see a 
clear difference between the two cases. While in the 10 case we need abont 10^ samples to reach 
a p-value close to 0.05, in the NO case in about one percent of the cases we have p > 0.05. The 
model needs a severe fine-tuning of the “anarchical” parameters in the 10 case, while the NO 
one is realized much more naturally. The most probable values of the profile parameters give 
Fg ~ (0.07,0.22,0.63), showing that approximate Anarchy is an outpnt rather than an input of 
the present construction. 

Focussing on the NO case, there is no preferred value of the leptonic Dirac CP phase. The 
lightest neutrino mass is predicted below 5 meV, corresponding to a hierarchical neutrino mass 
spectrum while lies in the range 0.1-5 meV, see fig. Any positive signal in the current 

generation of experiments aiming at measuring neutrino masses or |m^^| in the lab wonld 
essentially rnle ont the model. The hierarchy in the right handed neutrino spectrum is very 
pronounced and the corresponding mass distributions are peaked around 10® GeV, 10® GeV and 
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Figure 4. The predictions for the masses of RH neutrinos obtained for Xmin/^ ^ ™ 

of NO and tan/3 = 50, from ref. [56] . 


10^^ GeV, as shown in fig. 

In summary, fermion masses and mixing angles are well described by the map in eqs. ( |6|30|31D , 
in terms of input parameters of order one, the elements of the Y matrices. Such a map can be 
realized in several different frameworks and does not necessarily need an underlying symmetry. 
The setup is compatible with both SU(5) and SO(IO) grand unification and with the known 
solution to the gauge hierarchy problem. On the weak side, additional ingredients are probably 
needed to control the new sources of FCNC and CP-violations arising from new flavoured physics 
at the TeV scale. Moreover all entries of the Y matrices are independent free parameters and it 
is not possible to make absolute predictions, beyond the order-of-magnitude accuracy. This is 
clearly a major limitation, since we would like to test the theory at the level of the best available 
experimental precision. Finally the map in eqs. ( 6|30|31 ) might be an oversimplified version 
of a more accurate description and indeed there are several variants of the frameworks briefly 
mentioned in section 2 that modify the results of eqs. depoIsTI ). 


5. Flavour symmetries 

Given the successful use of symmetries in other contexts of particle physics, it is natural 
to investigate whether Yukawa couplings can be constrained by a flavour symmetry. The 
largest possible classical flavour symmetry of a theory with the particle content of the SM 
is Gmfv = U(3)^ and corresponds to the limit in which the Yukawa couplings are turned 
off. The observed fermion masses and mixing angles break Gmfv almost completely to a 
residual symmetry that includes the weak hypercharge and the combination B-L. Similarly, in 
any realistic model based on flavour symmetries, the flavour symmetry group Gf is broken. 
In predictive models the breaking is spontaneous and occurs through the VEVs of a set of 
scalar fields ip transforming non-trivially under Gf. The VEVs {(p) are either postulated or 
determined by minimizing a Gj-invariant energy functional V(ip). The Yukawa couplings become 
dynamical variables evaluated at the minimum of V{(p): y{{ip)/Af). A huge number of models 
can be constructed according to this set of rules, depending on the choice of Gj (global, local, 
continuous, discrete, abelian, non-abelian), and on the choice of representations for scalars and 
fermions. 

An attempt to start from the full Gmfv symmetry is described in refs. [571 Eg Eg Eol El]. 
The framework is that of MFV m- Under S'C/(3)^ quarks transform as in eq. (|24[). Yukawa 

















couplings are promoted to spurions transforming as 


Vu — (3, 3,1) Ud — (3,1,3) 


( 42 ) 


to ensnre invariance of the Yukawa interactions nnder Gmfv- By analyzing a general Gmfv- 
invariant potential depending on yu,d-, it has been proven the existence of stationary points 
where 

yu 4 oc diag(0, 0,1) Vckm = 1 , (43) 

pointing to an approximate U(2)^ symmetry of the qnark spectrnm and providing a good first- 
order approximation of quark masses and mixing angles. By extending Gmfv to the lepton 
sector, other stationary points have been identified: 


ye oc diag(0,0,l) 


rrip — UpMNS^’^^PMNS 


(44) 


rhjy = diag(mi, m, m) 


I 1 0 


UpMNS = 


0 

0 4 4 


V2 y/2 
f) —— -L 


X diag(-i,-i, 1) 


(45) 


Such a solution nicely exhibits maximal 023) bnt needs mi = m to further enforce a large, 
undetermined 012- Goldstone bosons arising from the breaking of Gmfv can be eliminated by 
gauging the flavour symmetry |62j . Correction terms are needed to promote the leading-order 
picture into a more realistic theory. If there is new physics close to the TeV scale, the advantage 
of assuming the largest possible flavonr symmetry Gmfv is the high degree of protection against 
FCNC induced by the new flavoured degrees of freedom. If the only sources of flavor symmetry 
breaking are the SM Yukawa couplings, such a maximal symmetry is minimally violated m- 
The classification of the expected effects is unambignons in the quark sector, and can be extended 
in several ways in the lepton sector [l6l |63] . 

Continuous flavour symmetry groups such as SO(3) and SU(3) have been proposed [6H[65l[^ . 
also in the context of grand unified theories |67l [68] , with the three fermion families assigned 
to an irredncible triplet representation. Charged fermions of the first two generations are much 
lighter than those of the third generation and consequently within SU(3)/SO(3) we are forced 
to introduce large breaking terms. Alternatively we can start from the smaller flavonr group 
U(2) and assign the first two generations to donblets and the third one to singlets. Within the 
simplest realisation, such an assignment in the quark sector leads to |69j 


yu,d — 


0 <4 0 \ 

^u,d ^u,d G{eii^d) 1 {yu,d)33 

0 Oi€u,d) 1 / 


(46) 


where phases have been omitted. To correctly reprodnce quark masses without appealing to 
cancellations among the matrix elements we need <C \iu,d\ 1) which corresponds to a 
sequential breaking of U(2): 

17(2) A 17(1) nothing . (47) 

The following relations can be derived [101 HI ] by a perturbative diagonalization of yu 4 '- 



( 49 ) 











up to corrections of relative order e, numerically close to few percent. The phase is arbitrary. 
These predictions can be translated into constraints on the (p, r/) plane and compared to the 
region presently allowed by the data, under the assumption that the SM correctly describes all 
relevant processes. Using as inputs quark masses, the maximally allowed values for sin 2/? and 
for derived from (49) are too small and the above set of Yukawa matrices are now 

excluded at the Scr level [zaizaiQ Modifications of the ansatz ( |46[ ), where either the 13 matrix 
element is non-vanishing or the element 32 is of the same order of the element 33, have been 
proposed to recover agreement with the data m- 

This is a nice example of a model of fermion masses based on symmetry requirements and 
leading to testable predictions. In this case the predictions were initially supported by data. 
Later on the precise measurements at the B-factories and the improved theoretical knowledge of 
both perturbative and non-perturbative QCD effects have considerably reduced the errors on the 
CKM matrix elements. Also the knowledge of light quark masses has improved and deviations 
from (49) are now signihcant. Unfortunately not all present models and ideas can be tested at 
the same level of accuracy. 

Independently of their dynamical origin, yu,d in eq. (46) are an example of textures, a simple- 
minded but predictive approach to the problem of fermion masses and mixing angles pioneered 
by Fritzsch [ZSIITT] and Weinberg [T|. In this approach we assume a privileged flavour basis 
and some special energy scale where some of the entries of the matrices yj vanish. In the 
quark sector, this can lead to relations between the elements of the CKM mixing matrix and 
the quark mass ratios which can be precisely tested. Not all zeros give rise to relations among 
physical quantities. By performing weak basis transformations, we can generate zeros in yj that 
have no physical implications [78]. For instance, working with two generations, we can always 
set to zero the ( 1 , 1 ) entry of both y^ and yu by performing a unitary transformation on the 
SU(2)-singlets quark fields. This transformation corresponds to a change of basis that has no 
physical consequences. If, in such a basis, we impose that yd and yu are symmetric matrices, this 
represents a non-trivial requirement. By further assuming small off-diagonal matrix elements, 
we get in first approximation the well-known Gatto-Sartori-Tonin relation m, eq. (Esj). Other 
unitary transformations, such as permutations, preserves the number of zeros and allow to divide 
the possible patterns of yd and into equivalence classes with identical predictions. It is always 
possible to enforce texture zeros in arbitrary entries of the fermion mass matrices by means of 
Abelian symmetries m- 

Special attention has been paid to hermitian textures, yu,d = Uu d- such texture is called 
n-zero if yu and yd have a total of n zeroes, the off-diagonal ones counting a half. A typical 
texture is 

/ _ 0 _ A/ 0 \ 

Vf = \ ] if = u,d) , (50) 

VO Bf Cf J 

where Du = D^ = 0 in the original proposal by Fritzsch. The 6 -zero Fritzsch texture is ruled 
out since it predicts a too large |I4fe|- More general hermitian 6 -zero and 5-zero textures have 
been analyzed under the assumption of hierarchy among the matrix elements . All possible 
combinations of hermitian 6 -zero and 5-zero textures are ruled out by now q because they 
cannot reproduce \Vch\, with the only exception of yj in (50) with Du 7 ^ 0 and D^ = 0, which 
is marginally compatible with present data [82|. The 4-zero texture in (50) with both Du and 


^ The most recent results for sin2/3 and |F„i,| (sin2/3 = 0.695 ± 0.023(0. 692to oi8), \Vub\ = 36.3± 1.2(35.7tl'.5) ^ 
lO”'^ [ 74 ] ([75])) are compatible with the ones in ref. [73] (sin2d = 0.739 ± 0.048, \Vub\ = (35.7 ± 3.1) x 10“"^) 
and, at the same time, more restrictive. 

^ According to ref. [5T], the hve independent hermitian 5-zero textures originally proposed by | 80| are still viable. 
However the values of d obtained from these textures in ref. m are too small, out of many sigmas from the 
currently allowed experimental range, except for the texture IV where the agreement is within about two sigmas. 









non vanishing is still allowed, provided the hierarchy in the 23 block is mild [831184| . In 
particular \Bf/Cf \ should be 0(0.1). In this case cancellations among the matrix elements are 
needed to get ms/rrib, rric/rnt and |I46| in the correct range. Notice that such cancellations 
were excluded by assumption in the case of of eq. (46), where the hierarchy between quark 
masses of second and third generations was attributed to the smallness of e. 

More general textures have been analyzed in ref. [85]. If no symmetry requirement is imposed, 
there are viable textures that correctly fit the data. However they do not show any predictive 
power with respect to any of the quark masses and mixing parameters. When are symmetric 
matrices, there are several independent 4-zero and 5-zero textures that allow to predict some of 
the light quark masses in terms of the mixing angles and the remaining masses. 

In the lepton sector, working in the flavour basis where the charged lepton mass matrix is 
diagonal, we can study textures of the symmetric matrix for light Majorana neutrinos | 86 | . Those 
with three or more zeros are experimentally excluded, while data still allow seven independent 2 - 
zero textures [871188] . The requirement that zeroes of the neutrino mass matrix should be present 
in the flavour basis can be relaxed and more general textures have been investigated. Both 
cases of Dirac and Majorana neutrinos have been analyzed and many independent textures in 
agreement with the existing data have been identified. For a systematic study of all possibilities 
of texture zeros in the charged-lepton and neutrino mass matrices and for a recent review on 
the subject we refer the reader to ref. |89| and ref. [90], respectively. 

Other choices of G/ have been considered. For instance refs. [9ll|92l|93| study the case of 
Gf = U{2)^, as an alternative to MFV to suppress FCNC in supersymmetric extensions of the 
SM. In that context, the advantage of Gf = U{2)^ is that squarks of the first two generations 
can be taken significantly heavier than the third generation ones. Furthermore the presence of 
separate U{2) for left and right-handed fields, provides a sufficient protection of flavour-violating 
effects in the right-handed sector. 


6. Discrete symmetries 

The data from neutrino oscillations before 2012 were supporting flavour symmetries, especially 
through the indication of a vanishing reactor angle 613 and a maximal atmospheric mixing 
angle 623 , features that are difficult to attribute to an underlying theory based on pure chance. 
Today we know with accuracy that ^13 is neither vanishing nor particularly small, its size being 
comparable to that of the Cabibbo angle. 

Recent global fits [Ml [Ml IM] (see table B favor a deviation of the atmospheric mixing 
angle from the maximal value by several degrees and, at the same time, show a preference 
for a maximal Dirac CP-violating phase Scp, though the whole range from 0 to 27r is still 
allowed at 3a. These two features are closely related, since they are mainly driven by the direct 
comparison between Pee measured by reactor experiments, which essentially determines ^ 13 , and 
the conversion probability P^e measured by MINOS and T2K, which in turn is sensitive to a 
combination of ^ 13 , 623 and 5cp. There is not yet a clear indication of the octant 623 falls in. 
Furthermore the most precise single experimental determination of the atmospheric angle HZ], 
6*23 = 0.514j'iQg5g(0.511j'iQg55) for NO (10), is still compatible with 7 r/ 4 . Probably it is too early 
to conclude that a maximal 623 is ruled out by data. 

There are few known mechanisms to generate a maximal atmospheric angle. We know that 
023 cannot be made naturally maximal by renormalization group evolution, barring a fine-tuning 
of the initial conditions or ad-hoc threshold effects [981 [99] irnniiToT] . Moreover, in the context 
of flavour symmetries, 623 is completely determined by breaking effects, if we accept that e and 
fL masses vanish in the limit of exact symmetry m- Indeed, if the masses of e, /x and r are 
all non-vanishing when the symmetry is exact, then large 0 ( 1 ) breaking terms are needed to 
explain the relative hierarchy among charged fermion masses. Thus it is more natural to assume 
that at least nie and are proportional to small symmetry breaking parameters. In this 







case, when the symmetry is exact, the lepton mixing matrix is determined up to an arbitrary 
rotation coming from the charged C/U lepton sector and the atmospheric mixing angle can only 
be determined when the symmetry breaking parameters are turned on. We have no examples 
of a maximal 622 , from the breaking of an abelian symmetry. If a nearly maximal atmospheric 
angle is not due to pure chance, we are left with broken non-abelian symmetries. 

Before the measurement of @13 a particularly attractive lepton mixing pattern was the 
tribimaximal one 
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Table 4. Results of global fits to neutrino oscillation data from ref. [Ml ESI [96] for normal 
[NO] and inverted [10] ordering. 
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The good agreement between TB mixing and pre-2012 data strongly supported the idea that 
the true mixing matrix could be described in terms of small corrections to a LO mixing matrix 
UpMNSi which could be derived from symmetry considerations. The simplest way to reproduce 
the TB mixing pattern is by exploiting discrete flavour symmetries |in3[ Iin4[ 11051 Iin 6 [ 1107] . 
The theory is invariant under a discrete flavour symmetry Gf, broken down in such a way that 
neutrino and charged lepton sectors have different residual symmetries Gu and Ge, at least in 
a LO approximation where small effects are neglected. If neutrinos are of Majorana type, the 
most general group leaving rriy invariant and the individual masses rm unconstrained is .^2 x -^ 2 , 
a finite group. The subgroup Ge can be continuous, but Ge discrete is the simplest option. 
We require a sufficiently large Ge to distinguish the three charged leptons. For instance we 
can choose Ge = Zn {n > 3) or Ge = Z 2 x Z 2 . Once Ge and Gy have been chosen inside 
Gf, the embedding automatically fixes the relative alignment of mjmi and niy in flavour space. 
Lepton masses are unconstrained but Upj^j^g is determined up to Majorana phases and up to 
permutations of rows and columns. This freedom apart, this setup predicts the three mixing 
angles and the Dirac phase S^p. In most concrete models, where symmetry breaking is 
achieved via VEVs of a set of flavons (p, the LO results are modified by small corrections of 
order u = {(p)/Af 

UpMNS = Upj^jpfg + 0{u) . (52) 

Before 2012 , in the specific case Upj^j^g = Utb these corrections were expected to be very 
small, of the order of few percent [IPS] 1109] . not to spoil the good agreement in the predicted 
value of the solar mixing angle. On this basis the simplest models reproducing Utb at the LO 
predicted ^13 not larger than few degrees, now proven to be wrong by experiments. Discrete 





















flavour symmetries can also be extended to quarks and even incorporated in GUTs, but in the 
existing constructions the symmetry has to be badly broken in the quark sector. Discrete flavour 
symmetries are also relevant in the so called indirect models mi- In this case the breaking 
oi Gf leaves no residual symmetries and the role of the flavour group is mainly to get specific 
vacuum alignments of the scalar fields that control fermion masses. 

Several modifications of the simplest models based on discrete symmetries have been proposed 
to match the most recent data. If we keep adopting = Utb as LO approximation, 

an economic way to reproduce the actual value of @13 is to introduce large correction terms, 
0{u) Rs 0 . 2 . This is also viable in other schemes where differs substantially from Utb, 

such as the so-called bimaximal mixing. 
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(53) 


Introducing large corrections has the disadvantage that beyond the LO the number of 
independent contributions to the mixing matrix is generally quite large. If their typical size 
is about 0 . 2 , all mixing angles tend to be affected by generic corrections of this type and 
predictability is lost. Moreover large correction terms are dangerous if new sources of flavour 
changing and/or CP violation are present at the TeV scale. Thus some assumptions about the 
dominant source of corrections are needed. For example, a reasonable possibility is to perturb the 
BM mixing pattern by a rotation U 12 from the left, possibly originating from the diagonalization 
of the charged lepton sector |110l lllll I112j 


UpMNS = Ui2{a,S)UBM 
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To first order in a we have 
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(55) 


By eliminating (a,<5) we get a relation between sin^ 0 i 2 , sin 0 i 3 , and bcp, plotted in fig. This 
model predicts dcp close to tt in order to reproduce correctly sin^ 0 i 2 , as can be seen from eqs. 

Another possibility is to relax the symmetry requirements. S 4 is the smallest group 
reproducing TB mixing through the breaking down to = Z 3 and G^ = Z 2 x Z 2 , 
whose generators are T and {S,U), respectively]^ In the basis where T and the charged 
leptons are diagonal, the element U coincides with the so-called exchange symmetry 
[ml [112 Ha [Ml ttni [mi, directly responsible for the vanishing of 0 i 3 and for 023 being 
maximal. We can avoid having 0 i 3 = 0 and 623 = vr/d if is a single Z 2 subgroup generated 

^ It is well-known that, in concrete models, the TB mixing pattern can also be obtained from the group A 4 , 
generated by S and T, if [/ arises accidentally due to a particular field content. 






















Figure 5. Left panel: contours of equal sin^ 0 i 2 in the plane (sin 6 * 13 , Jcp/^)) when Upmns = 
U12 X Ubm- The brown region is excluded at 3cr, assuming [0.0177 0.0294], [0 27 r ] and 

[0.278-^0.375] as 3a ranges for sin^ia, 6 cp and sin^ 0 i 2 , respectively. Right panel: contours of 
equal sin^ 6*23 in the plane (sin 0 i 3 , 6 cp/tt), when U pm ns = Up M2 ■ The brown region is excluded 
at 3c7, assuming [0.0177 0.0294], [0-^ 27r] and [0.392-^0.643] as 3a ranges for sin 0 i 3 , 5cp and 

sin^ 023 , respectively. 


either by the element S or by the element SU. When the preserved parity is S, the mixing 
pattern, TM 2 , is trimaximal and corresponds to 

( cos a 0 sin a \ 

0 10, (56) 

sin a 0 cos a / 

with both a and 6 unconstrained. When the preserved parity is SU, the mixing pattern, 
TMi, is also of trimaximal type and is given by UpMi = UpB U 23 {a, 6 ), where U23{a,5) is the 
transformation analogous to Ui 3 {a, 6 ), acting in the 23 plane. The mixing angles and the Dirac 
phase are predicted in terms of {a, 5) and we get two relations among physical quantities, shown 
in Table |119l 11201112111122111231112411125L I126j . The interesting feature of these relations is 
that the predicted deviations from TB are linear in sin 0 i 3 for sin^ 623, and quadratic for sin^ 0 i 2 , 
known with much better precision. One of these relations is plotted in fig. [^in the case of TM 2 , 
from which we see that a substantial improvement in the data is needed to test this possibility. 
Explicit models based on ^4 realizing the TM 2 breaking pattern were indeed proposed before 


Table 5. Sum rules for TMi ^2 mixing patterns. 


TMi TM2 

sin^ 012 = 5 - I sin^ 0i3 -|- 0(sin^ 013) sin^ 0i2 = g + 5 sin^ 0i3 -|- ©(sin"^ 013) 

sin^ 023 = 5 - \/2 sin 0i3 cos 5 cp + 0(sin^ 013) sin^ 023 = 5 + sin 0i3 cos 5 cp + 0(sin^ 013) 


the measurement of 0 i 3 |127| . The possibility of reducing the residual symmetry Gy to Z 2 can 
be systematically investigated mi- 







































A further possibility is to look for alternative LO approximations where @13 is closer to the 
measured value. Remarkably, an infinite set of groups Gf giving rise to LO approximations 
closer to the data has been found. Of particular interest is the special form of trimaximal 
mixing TM 2 in (56), with both a and 6 quantized, reproduced by groups of the series A( 6 n^) 
|129l 11301113111132j . For example, choosing n = (4, 8 ,10) we have a = (± 7 r/ 12 , ±7r/24, ±7r/15) 
and sin^^j’g = (0.045,0.011,0.029). The Dirac phase is zero (modulo vr). In fig. § the values of 
I f 7 e 3 1 are plotted versus n |132j . Other discrete groups have been investigated in ref. |133j . Very 


V13I 



Figure 6. Possible values of \Ue 3 \, indicated as |Vi 3 | in the vertical axis, versus n in A( 6 n^) 
models, from ref. m- The lines denote the present approximate 3a range of \Ue 3 \- Examples 
include |Vi 3 | = 0.211,0.170,0.160,0.154 for n = 4,10,16,22, respectively. Each value of Ue 3 
allows for two values of 623 with 5cp = ^ and 6 cp = p given by 623 = 45° ^di 3 /y /2 respectively. 

remarkably, a complete classification of all possible mixing matrices \Upj^j^g\ generated from 
any finite group has been recently carried out in ref. m- 

Another development consists in combining discrete and CP symmetries |135( I136j and 
exploring the symmetry breaking patterns such a combination can give rise to. A well-known 
example is that of the so-called fir reflection symmetry jl37l 11381113911140] (not to be confused 
with the exchange symmetry), which exchanges a muon (tau) neutrino with a tau (muon) 
antineutrino in the charged lepton mass basis. If such a symmetry is imposed, the atmospheric 
mixing angle is predicted to be maximal, while ^13 is in general non-vanishing for a maximal 
Dirac phase 6 and the Majorana phases vanish. The solar mixing angle remains unconstrained. 

A general formalism which combines CP and flavour symmetries m can be used to 
constrain the lepton mixing matrix. A theory symmetric under CP and under a discrete 
flavour group Gj is assumed to have residual symmetries G^, generated by some elements Qi 
and Cl, = Z 2 X CP, generated by a parity Z and a CP transformation X. The action of 
X in flavour space can be non-trivial m and should respect a set of consistency conditions 
[135111421114311144] . The residual symmetries Gg and G^, imply the following conditions on mjmi 
and m^: 

Q\{m\mi)Qi = {m\mi) , Z'^m^Z = mi, , Xm^X = ml . (57) 

These conditions are strong enough to determine Upj^j^g completely, up to one real parameter 
9, ranging from 0 to vr: 

^PMNS — ^PMNsiQii 0 < 0 < vr . (58) 

Mixing angles and phases, both Dirac and Majorana, are then predicted as a function of 9, 





























































































Figure 7. Results for the mixing parameters sin^is, sin^ 0 i 2 for Case I (straight line) and 
Case IV (dashed line), from ref. [135] , We mark the value 9-^^ of the parameter 6 for which the 
functions have a global minimum with a red dot. 3 ct ranges for the mixing angles are also 
shown 


modulo the ambiguity related to the freedom of permuting rows and columns and to the intrinsic 
parity of neutrinos. The formalism is completely invariant under any change of basis in field 
space. The physical results only depend on the initial symmetry and the residual symmetries 
specified by (Qj, Z^X). An interesting example is provided by Gf = ^ 4 . An exhaustive analysis 
has been presented in ref. |135| . The residual symmetries can be chosen as Ge = Z^, generated 
by the element T, and G^ = Z 2 x GP, generated by {Z,X). The parity transformation Z can 
be either S (case I) or SU (case IV) and a consistent GP transformation X acting on the lepton 
doublets coincides with the /i — r reflection symmetry in the basis where T (and the combination 
mjmi) is diagonal. Thus the predicted mixing pattern has a maximal atmospheric mixing angle, 
a maximal Dirac phase, vanishing Majorana phases and there is a relation between the solar 
angle and the reactor angle, shown in fig. 

Recently several explicit models combining CP and flavour symmetries have been proposed 
and several series of discrete groups have been investigated in combination with CP |145l I146L 
1147111481114911150j . Other approaches making use of CP and discrete symmetries are described 

in [mlfT5mfT5mfTMfT55l[TM[TFn[T^[TM[TMfT6TlfT62lfT63] . 


7. Hints from empirical relations 

Empirical relations among fermion masses and/or mixing angles have been frequently suggested 
as a clue towards a solution of the flavour puzzle. Here as an example I will comment one of the 
most striking ones, Koide’s formula for charged lepton masses HMj: 



y/m-e + m^ + mr 


( 59 ) 









































The experimental values of the pole masses rUe, rrifj,, rur are mi¬ 


me = 0.510998928(11) MeV 
= 105.6583715(35) MeV 

mr = 1776.82(16) MeV . (60) 


By using as input me and in eq. (60), from Koide formula we find m,- = 1776.97 MeV, with no 
appreciable dependence on the experimental errors of me and m^ and in perfect agreement with 
the measured value. There are several attractive features of Koide formula, beyond its simplicity. 
It is independent of the choice of mass units. It can be expressed as K{me, m^, mr) = 0, where 
K is an homogeneous function of the charged lepton masses, symmetric under permutations of 
me, m^, mr- This makes the formula particularly attractive and has stimulated lot of interest 
and of activity aimed at deriving or embedding it within a more fundamental framework. An 
unsatisfactory feature of the formula is that it requires as inputs pole masses and not running 
masses, which would be more adequate if masses are believed to originate from a common scale 
A f. This can be seen in table where the running Yukawa parameters of the charged leptons 


Table 6. Values of the running SM Yukawa couplings for the charged fermions in the MS 
scheme, at selected renormalisation scales /r, from ref. jl66j . The fourth and the fifth rows 
show the prediction of the mass relations in eqs. (59) and (62), respectively, using as inputs the 
first two rows, yeifJ-) and The errors indicated in brackets affect the last digit and are 

dominated by the uncertainty in y^{n). 



// = Mz 

y = 1 TeV 

y = 3 TeV 

fi = 10 TeV 

ye / 10-® 

2.794745;n[!om6 

0 Q/1 qoH“ 0.0022 

^.545^q 0021 

2.8646t[J:“i 

9 q'7CoH“ 0-0042 
z.o/oz_0 0039 

2//. / 10-^ 

5.899863;noom8 

6 0127+°-°°^^ 

QQ44 

6 . 0473;°:“67 

0-U;Di_Q Qgg2 

yr / 10-2 

1-OU290U_Q QQQQQ^ 

1 0991 

i.uzzio_o_oo077 

1.0280 ±0.0011 

i.uozy_g 0045 

2/f/10-2 

0.990448(3) 

1.0094(7) 

1.015(1) 

1.020(1) 

yP / 10-2 

0.991610(3) 

1.0106(8) 

1.016(1) 

1.021(1) 


in the MS scheme are listed for several values of the scale /r. The fourth row shows the value 
of the r Yukawa coupling y^{y) as derived from the Koide formula using as inputs ye{y) and 
yfiin)- In the range /r = 0.1 10 TeV the mismatch between y^{y) and yr{y) is close to one 

percent, while the accuracy with which yr{y) is estimated in ref. |166j ranges from 10“'^ to 
10“^, thus making the discrepancy significant. For example the value of the tau mass at the 
scale Mz predicted by the Koide relation is smaller by about 20 MeV, while the experimental 
precision on the same parameter is close to 0.2 MeV. This gap mainly arises as the effect of the 
pure QED running going from me to mr as can be seen from the leading order relation between 
pole masses and MS masses: 


mi{fi) = mi 



(l + ^log 



(61) 


At energies higher than 1 10 TeV, the success of the Koide relation depends on unknown 

physics. If the SM merges into a supersymmetric theory, the evolution of the Yukawa couplings 























depends on additional parameters, like the superparticle thresholds and tan/3. Just above the 
superpartners mass threshold the Yukawa coupling are known to a much worst precision, around 
one percent m, and the Koide relation might again be compatible with the extrapolated data. 
At the GUT scale, larger uncertainties are induced on ye,/i,T(TJG[/T) by tan/3. 

More in general, we can ask what is the probability of finding a simple homogeneous relation 
among charged fermion masses holding to an accuracy similar to the one of the Koide relation 
at a given scale /r. It is difficult to answer quantitatively this question, but I personally think 
that such probability is relatively high. As an example consider the following relation, valid for 
charged fermion masses rather than for their squared roots: 

unie + + nir 

n^e + nT'fj. + W-r 

This relation, symmetric under permutations of the flavour labels, produces the outputs in the 
fifth row of table We see that the predictions of are as good as the one from the 

Koide relation. Probably it is not surprising that scanning thousands of possibilities we can find 
“simple” relations working at the level of the percent accuracy. Much more difficult is to find, 
through these relations, a direct link to some unknown fundamental layer of particle physics. 

Conclusion 

We are witnessing a continuous experimental progress in flavour physics. In neutrino physics, 
squared-mass differences and mixing angles are known to an accuracy that approaches the 
percent level. The reactor angle ^13 is away from zero by many standard deviations. For the 
first time global fits hint at a non-trivial Dirac phase. While the new data have been effective in 
ruling out many models of fermion masses and mixing angles, as a matter of fact no compelling 
and unique theoretical picture has emerged so far. 

Present data can still be described within widely different frameworks. Based on our 
experience with gauge interactions we might hope that the flavour sector becomes simple and 
symmetric at a high energy scale, with a small number of relevant parameters providing a 
complete description. It is fair to say that we have not been able to identify a clear symmetry 
pattern from data so far. Before 2012 discrete symmetries were considered as a good candidate. 
In particular those predicting a nearly tri-bimaximal mixing were favored by data, but the 
prediction of @13 turned out to be wrong. The evidence for discrete symmetries in the quark 
sector is very poor and in a unified description of all fermions this kind of symmetry is typically 
badly broken in the quark sector. The whole approach is too much centered on the lepton 
mixing properties, while a description of the fermion masses seems to need additional ingredients. 
Several modifications of the simplest schemes to accommodate the present data are still possible 
and have the advantage of being quantitatively testable. But the real open question is whether a 
non-trivial implementation of discrete symmetries exists encompassing quark and lepton sectors 
in a unified picture and providing a quantitative description of both masses and mixing angles. 

There are other models where fermion masses and mixing angles are mapped into a large 
number of irreducible and unconstrained order-one parameters, thus incarnating the Anarchy 
idea. For their intrinsic nature these models essentially escape experimental tests going beyond 
the order-of-magnitude accuracy. However we cannot fail to be impressed by the fact that they 
can provide a common description to both fermion masses and mixing angles, that they are 
compatible with grand unified theories and that they can be derived within widely different 
theoretical frameworks. The fact that this kind of models can be implemented even in a highly 
constrained setup such as an SO(IO) grand unified theory is really remarkable. As a drawback, 
in these models the bounds on the scale of new flavored physics is typically pushed above the 
10 TeV scale, reducing the possibility of testing these ideas. 
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(62) 






Flavour remains a fascinating mystery, still eluding all our attempts to find the rationale 
underlying our observations. Has this puzzle any solution? Are we misled by the questions 
we have formulated so far? Will we ever have access to the flavour scale? Man has long been 
fascinated by the mystery of planet motion. Surprisingly precise measurements and accurate 
predictions already existed in remote epochs. For a long time the most reliable models were 
based on the special character of geometrical figures like the sphere or the circle. Attempts 
to explain the relative sizes of the solar planetary orbits revealed themselves misleading. More 
accurate observations, perseverance in identifying the correct questions and renunciation of old 
prejudices opened new perspectives to the scientific thought. Will this be the fate of the flavour 
puzzle too? 
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